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While numerous effective methodologies exist for displaying the magnitude 



of individual vectors, the ability to observe commonality of individual vector orientation 



within vector fields containing thousands of individual vectors on a surface or within a 



volume can be difficult to discern. 



5 



The applicant is aware of several attempts in prior art, which relate 



particularly to the visualisation of three-dimensional vector fields. Reference may be had by 



turning to the following publications h e reinaft e r furth e r r e f e rr e d to in 



sup e rscript 



tForm 



PTO/SB/08B where references are cited thereon. 

Reference 4 ': "Visualizing th e Magnitud e and Direction of Three - Dim e nsional 
10 V e ctor Fi e lds'' by Timothy J. Peter in AP - S International Symposium Digest, IEEE 
Ant e nnas and Propagation Soci e ty, Volum e 3, 1993, pag e s 1 1 18 to 1121. 

R e f e r e nc e'- — "Visnnlizntion — b£ — Thr ee- Dim e nsionnl — V e ctor — Fields — m4 
Functions" by J. Simkin, IEEE on Visual 3D Fields Digest, 1995/205, pages 1/1 2; and 

Reference*: "V e ctor Fi e ld Animation with Maps" by Yamrom, IEEE 
15 Computer Graphics and Applications, Volum e 15, Part 2, pag e s 22 to 2 4 , March, 1995; 

R e f e r e nce 4 : "Combin e d 3D Int e rpretation of Airborn e , Surfac e and Bor e hol e 
Vector Magnetics" of Mueller, E. L, Morris, W. A., Killccn, P.G., and Balch, S. at the 
McConncll Nick e l Deposit, Proc ee dings of Exploration 97: Fourth D e c e nnial International 
Conf e rence on Mineral Exploration, e dit e d by A. G. Gubins, 1997, p. 657 - 666. 
20 R e ference*: "Mapping R e mnant R emanent Magn e tism Using the Local Phase" 

of Thruston, J. B., in Geophysics, Volume 66, No. 1, July August, 2001, p. 1082 1089. / 

Current methodologies for indicating vector orientation consist of 
streamlines, hedgehogs, stream-ribbons, stream-polygons and textures maps 3 . In all these 
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methods there are difficulties in observing the vector orientation orthogonal to the 
observation surface or in selection of unique desired orientations. Figure 5 of the drawing 
sheets illustrates a vector grid that shows this effect. A novel texture mapping procedure is 
hereby disclosed that overcomes many of the shortfalls of the earlier methods. 
5 Many applications exist for the use of an effective display of vector 

orientation including studies in fluidics, electromagnetic fields, finite element analysis, 
magneto hydrodynamics, plasma physics, astrophysics magnetic and digital elevation 
models (DEMs). We will present an embodiment in which the application of the present 
method is illustrated in the field of geophysical geomagnetic mineral exploration. It is to be 

10 recognized that this embodiment is but only one example of the usefulness of this Vector 
Orientation Visualization Method. 

Geomagnetic surveys have many uses and are of particular interest in mineral 
exploration, as they may identify areas likely to contain some types of valuable deposits. 
They can also identify various lithologies, and provide some types of structural information. 

15 A common method of performing such a survey is to measure the magnetic 

field at various points along the Earth's surface while recording the magnitude variations in 
the magnetic field strength. Although dip needle surveys have been used in the past, until 
recently, the orientation of the magnetic field was not generally measured. The magnitude 
of variations in the readings was primarily due to variations in the magnetic susceptibility of 

20 the constituent minerals within the rock, and to a much lesser degree, the variation of the 
orientation of the r e mnant remanent field. Such surveys, while giving useful data, which 
related primarily to lithological effects, disregarded the information of the orientation of the 
magnetic field. 
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Recently, due to the introduction of tri-axial gradient magnetometer surveys, 
it has become possible to retrieve vector oriented gradient data. In addition, tri-axial 
fluxgate total magnetometer surveys can also be used to measure the orientation of the total 
magnetic field vector. Tri-axial fluxgate magnetometers have also been used in along a 
5 borehole to effectively delineate the partial outline of a dipole about a ferromagnetic ore 
body 4 . 

The ability to effectively display the orientation of the magnetic vectors over 
a surface or within a block of rock as by borehole surveys can greatly add to the geologic 
understanding of the survey area and can lead to a better understanding of the location and 

10 extent of several varieties of ore deposits. 

To understand the benefit of knowing the orientation of the vectors within the 
survey area it is important to understand the means by which variations in the orientation of 
the vectors can occur. After a high temperature thermal event, certain minerals in the earth 
will heat enough beyond their Currie Point Ne'el temperature to lose their magnetism. 

15 When the temperature lowers, these minerals regain their magnetism and will now have a 
thermal femnatt tremanent magnetic field aligned with that of the earth at the time of going 
through their Currie temperature. Since the earth's magnetic field is known to change 
directions and flip at various times, the differently oriented r e mnant remanent magnetic 
vectors of these remagnetized minerals can indicate spatially segregated intrusives of 

20 contemporaneous age, multiple pulsed intrusives within a similar rock type, the extent of 
their thermal aureoles and the presence of folds and faults. Similarly alteration due to 
chemical -known as chemical remttant remanent magnetism, CRM- and pressure changes - 
known as piezoreroremi remanent magnetism- can result in the formation of new minerals, 
which will retain a femnantremanent orientation of the earth's field orientation at that time - 
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known as detrital femnani remanent magnetism (DRM) caused by the quiescent 
sedimentation of fine magnetic material aligning with the earth's magnetic field at the time 
of sedimentation.- DRM is another remnance whose variation can be used to recognize later 
faults and folds. 

5 For the purpose of this embodiment a magnetic survey can be performed 

using either a tri-axial fluxgate magnetometer or tri-axial gradient magnetometer. The 
vertical gradient can also be calculated if both horizontal gradients and the total magnetic 
field are known. These types of surveys produce the magnitude of the magnetic field for 
each of the three-orthogonal directions. The data is collected for process and analysis over 

10 the survey area by establishing a two-dimensional grid over the area of interest. 
Measurements of the magnetic field are taken at specific points on the grid either by aircraft 
or by an operator on the ground working along a cut or demarked grid. . Each collection 
point is identified as a point on the two-dimensional grid, and will have an associated three- 
dimensional vector representing the magnetic field at that point. Alternatively, if the data 

15 should be collected along a borehole as a function of its distance from the collar, the location 
of the three-dimensional vector can be represented from the surveyed borehole information 
as a point in 3D space with x, y, z co-ordinates. The x,y co-ordinates can match the two 
dimensional grid of the surveyed area of interest mentioned earlier. 

Peters Reference 1 : This reference describes a method of representing an 

20 electromagnetic vector field defined over a three-dimensional surface. A triangular grid 
represents the surface, and the method uses color shaded magnitude contours and directed 
vector lines to characterize the vector field. Referring to the Figures, it may be seen that in 
this reference, the direction of a vector field is represented using arrows, and the magnitude 
of the vector field is displayed using colors. Such visualization is suited to vector fields for 
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Therefore Xcaiibrated equals X m - Xbasestation(time) which is equal to X m - |Xbavg - Xb(time)|. Values 
for Ycaiibrated and Zcaiibrated are calculated in a similar manner. 

After the values Xcaiibrated, Ycaiibrated and Zcaiibrated for each measuring station 
have been calculated. A second step may be performed to remove the geomagnetic field 
5 vector from these values in order to obtain a residual magnetic value from the magnetic field 
data. 

This is achieved by calculating the Earth's magnetic field using the standard 
International Geomagnetic Reference Field (IGRF) procedures to determine the theoretical 
values with respect to the area of interest at the time of the survey These theoretical values 
10 obtained may be represented as Xigrf, Yigrf and Zigrf. 

The residual magnetic data values Xresiduai, Yresiduai and Zresiduai may then be 

Calculated from the equations Xcaiibrated- Xjgrf, Ycaiibrated - Yigrf, and Zcaiibrated - Zigrf, 

respectively. 

After calculating the residual magnetic values, further methods known to the 
15 trade 5 may be performed which causes the induced magnetic field value to be removed from 
the residual magnetic data values to obtain r e mnant remanent field values at the measuring 
stations. 

In order to calculate the fefflnan tremanent values X^^remm^u Y f ^ mft ««tremanent 

and Zfemi*Htf remanent, the equations Xresiduai - * Xindueed, Y re sidual — Y m duced, and Z res idual ~ Zindueed, are 

20 respectively determined. The fefflnan tremanent value represent the magnetic field data 
values were frequently acquired at the time the ferromagnetic rock forming minerals were 
formed or went through a Currie temperature episode of its ferromagnetic minerals. This 
femffiffi tremanent vector data can therefore allow for discrimination. After calculating either 
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the residual or r e mnant remanent data these x, y, z Cartesian values are transformed to 
mathematical spherical co-ordinates using the following equations: 

rmath = squareroot(X 2 + Y 2 + Z 2 ) where r is the radius of the unit sphere used 
to determine spherical co-ordinates; 
5 0(theta)math = arctan (Y/ X) where theta represents the azimuth angle, which 

is measured from the positive x-axis toward the positive y-axis theta and is a value between 
0° and 360°; 

<D(phi)math = arctan ((squareroot(X 2 + Y 2 )/ Z) where phi is the angle measured 
from the vertical positive z axis counterclockwise toward the x,y plane and is a value 
10 between 0°to 180°. 

These mathematical spherical co-ordinates are illustrated clearly in Figure 4. 

The mathematical spherical co-ordinates of the magnetic vector at each 
measuring station may then be represented by (rmath, ©math. ^math). 

A conversion of the mathematical spherical co-ordinates Ow™, ©math, Q>m&\h) is 
15 then translated into the geological co-ordinates (r g , @ g , Og) by the following equation: 
(r g , @ g , Og) = [(rmath), (©math - 270), (90- Omath)] as clearly illustrated in Figure 4a. These 
values are tabled in Figure 8. 

After calculating the geological spherical co-ordinates for each measuring 
station, a unique colour based legend may be determined by assigning colour or patterns to 
20 the spherical co-ordinates. Several colour models or notations exist which have direct 
symmetry with respect to spherical co-ordinates such as red, green, blue (RGB); 
Commission Internationale de l'Eclairage (CIE), hue, saturation, value (HSV); hue, 
saturation, luminance (HSL), CIE XYZ, YIQ, Munsell, TekHVC, CIE LUV etc... The last 
three of these models allow for a more continuous perception of colour to the human eye. 
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Additionally the magnitude of the vectors may also be shown on the coloured 
illustration either as dashed black and white line contours or as a Digital Elevation Model. 

Alternatively to the ground acquisition method described earlier, an aircraft 
flying over the grid may retrieve the above information. After the grid has been determined, 
5 the aircraft flies over a prescribed calculated grid lines to record the magnetic vector 
information at regular time intervals corresponding to a specified distance given the airspeed 
flown. In order to ensure that the aircraft follows the calculated grid lines, a GPS is be-to be 
used to guide the aircraft. In this aerial data collection no baseline or ground demarcation of 
lines is required. 

10 The magnetic field vector information is collected via instruments on the 

aircraft as they fly over the measuring station and stored on record keeping equipment, such 
as a data logger, located on the aircraft. Calibration and corrections to the data are carried 
out later after each day's flying. 

An example of the 3D versatility of this colour coded vector orientation 

15 methodology is also given in Figure 12 where the method illustrates vector orientation 
sample points collected along a borehole within a given volume by the use of coloured 
voxels. 

20 



